Testing Fundamental Physics with Distant Star Clusters: 
Analysis of Observational Data on Palomar 14* 



K. Jordii'2, E. K. GrebeP, M. HiW, H. Baumgardt^ M. Frank^, P. Kroupa^ H. Raghi^'^, 

P. Cote^, and S. G. Djorgovski^ 



ABSTRACT 

We use the distant outer halo globular cluster Palomar 14 as a test case for classical vs. 
modified Newtonian dynamics (MOND). Previous theoretical calculations have shown that the 
line-of-sight velocity dispersion predicted by these theories can diS^er by up to a factor of three 
for such sparse, remote clusters like Pal 14. We determine the line-of-sight velocity dispersion 
of Palomar 14 by measuring radial velocities of 17 red giant cluster members obtained using 
the Very Large Telescope (VLT) and Keck telescope. The systemic velocity of Palomar 14 is 
(72.28 ± 0.12) km s'^ The derived velocity dispersion of (0.38 ±0.12) km s'^ of the 16 definite 
member star s is in agreement with the theoretical prediction for the classical Newtonian case 
according to iBaumgardt et al.l (|2005l ). In order to exclude the possibility that a peculiar mass 
function might have infiuenced our measurements, we derived the cluster's main sequence mass 
function down to 0.53 Mq using archival images obtained with the Hubble Space Telescope. We 
found a mass function slope of a = 1.27±0.44, which is, compared to the canonical mass function, 
a significantly shallower slope. The derived lower limit on the cluster's mass is higher than the 
theoretically predicted mass in case of MOND. Our data are consistent with a central density of 
po = 0.1 Mqpc~^. We need no dark matter in Palomar 14. If the cluster is on a circular orbit, 
our spectroscopic and photometric results argue against MOND, unless this cluster experienced 
significant mass loss. 

Subject headings: globular clusters: individual (Pal 14) - gravitation - stellar dynamics 
1. Introduction 
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Is classical Newtonian dynamics valid on all 
scales? On Earth classical Newtonian dynamics 
describes all non-relativistic phenomena very well. 
With the exploration and study of the Universe, 
we can explore low acceleration regimes that can 
not be studied in our backyard and we observe 
deviations from the expected classical Newtonian 
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behavior. E.g. the rotation curves of spiral galax- 
ies do not show the classically expected Keple- 
rian fall-off, b ut stay flat in the ou ter parts of 
these galaxies ( Sofue fc RubinllioOl ). These flat 
rotation curves are commonly explained by intro- 
ducing dark matter (DM). In the outer parts of 
the galaxies, DM is more abundant than regu- 
lar baryonic matter and the gravitational effect of 
the D M on the baryons results in a flat rotation 



curve (jRubin et al.l 119821 ). A major problem DM 



theory has encountered recently is the discovery 
that young tidal-dwarf galaxies also have rotation 
curves that imply a significant invisible matter 
component although they cannot be dominated by 
non-baryonic DM suggesting a n on-classical phys- 
ical solution ( Gentile et al ] l2007l) . 

An alternative theory to D M is modi f ied Ne w- 
tonian dynamics (MOND; [Milgrom I Il983al lbl: 



Bekenstein fc Milgromi Il984l) . In MOND, the 



flat rotation curves of galaxies can be fitted with- 
out any assumption of unseen matter. According 
to MOND, Newtonian dynamics breaks down for 
accelerations lower t han an ^ 1 x 10^ ^ cm s^^ 
((Begeman et a l. 1991: iSanders fc McGaug h 2002). 
The acceleration a in MONDian dynamics is given 
by the (heuristic) equation: 



ao 



a = On, 



(1) 



where ^J.{x) is an arbitrary function with the fol- 
lowing limits: 



if a; < 1 
if a; > 1. 



(2) 



Here, oat is the standard Newtonian accelera- 
tion and ao is the acceleration limit below which 
MOND is applicable. 

It has been claimed that MOND has difficulties 
explaining the merging of galaxy clusters, where 
the baryonic matter is clearly separated from the 
gravitational mass, as fo und by gravi tational lens 
l2006l) 



However, lAngus et al 



ing (Clowe e t al 

([20061 120071) demonstrated that such systems are 
consistent with MOND, but do require the exis- 
tence of some hot dark matter. 

Baumgardt et al. ( 20051 BGK05) proposed to 
use distant, outer halo globular clusters (GCs) to 
distinguish between classical and modified Newto- 
nian dynamics. They calculated the line-of-sight 



velocity dispersion for 8 Galactic GCs in classical 
and in modified Newtonian dynamics. For these 
GCs the external acceleration due to the Milky 
Way and the internal acceleration due to the stars 
themselves are significantly below the critical limit 
of ao- The expected velocity dispersions in case 
of MOND exceed those expected in the classical 
Newtonian case by up to a factor of three (see Ta- 
ble 1 in BGK05). 

Palomar 14 (Pal 14) is a diffuse GC in the 
outer halo of our Galaxy. Pal 14's sparseness, 
faintness, and large distance to the Sun makes 
it a difficult observational target, and therefore 
it did not receive much attention. The first ra- 
dial velo city for a Pal 14 member s tar was mea- 
sured by I Hart wick fc Sargent (Il978l) resu lting in 
81 ± 3 km s-^ lArmandroff et all |l992[ ). based 
on radial velocity measurements for two stars, re- 
ported a systemic velocity of 72 ± 3 km s~^. 
The deepest ground-based color-mag nitude dia- 
gram (CMD) of Pal 14 was published bv lSaraiedinil 
(|l997t) . He concluded that Pal 14 is 3-4 Gyr 
younger than halo GCs with a similar metallicity. 
Hilkcr (2006, H06) published photometric data 
on three GCs from the BGK05 sample: AM 1, 
Pal 3, and Pal 14. H06 confirmed Pal 14's youth 
of ~ 10 Gyr. The data from his study are used 
here to obt ain targets f o r our spectroscopic ob- 
servations. iDotter et al. (2008) published a pho- 
tometric study of Pal 14 based on archival data 
obtained with the Wide Field Planetary Camera 
2 on board the Hubble Space Telescope. The au- 
thors confirm Pal 14's relative youth. Here the 
same data are used to obtain the cluster's mass 
function. 

This paper is the second in a series that inves- 
tigates theoretically and observationally the dy- 
namics of dis tant star clusters. In the first paper 
(|Haghi et al l, 2009. HBK09), we derived theoreti- 
cal models for pressure-supported stellar systems 
in general and made predictions for the outer-halo 
globular cluster Pal 14. In the current paper, we 
present a spectroscopic and photometric study of 
Pal 14, as a test case for the validity of MOND. 
I.e. we are measuring the velocity dispersion of 
Pal 14 in order to compare the measured value to 
the predicted values made for MOND and classical 
dynamics by HBK09. Further, we are determin- 
ing the mass function of Pal 14 in order to infer 
the cluster's mass. The derived mass and velocity 
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dispersion are then compared to the predictions 
made by HBK09 for Pal 14 on a circular orbit in 
MOND. 

The paper is organized as follows: In Section [5] 
we describe the observational material. In Sec- 
tions [3] & |4] we present stellar radial velocities, the 
color-magnitude diagram and the mass function of 
Pal 14. In Section [5] we discuss the effects of our 
result for MOND and classical Newtonian grav- 
ity. The last Section concludes the paper with a 
summary. 



2. Observations and Data Reduction 

To distinguish between MOND and classical 
Newtonian dynamics we used two different kinds 
of observations. In order to measure Pal 14's 
velocity dispersion we obtained high-resolution 
spectra of red giant candidates towards Pal 14 
with the Ultraviolet-Visual Echelle Spectrograph 
(UVES; iDekker et all I2OOOI ) at the Very Large 
Telescope (VLT) of the European Southern Obser- 
vatory (ESO) in Chile and with the High Resolu- 
tion Echelle Spectrograph on the Keck I telescope. 
To be able to measure the cluster's mass function 
we used imaging data from the Hubble Space Tele- 
scope archive. In the following subsections we de- 
scribe the reduction process of our observational 
data. 

2.1. Spectroscopy Vifith UVES 

The photometry published by HOG shows the 
red giant branch and horizontal branch of Pal 14. 
Based on this photometry, we selected 16 of the 
17 brightest red giants of Pal 14 for spectroscopy 
with UVES at the VLT. Our target stars cover 
the magnitude range V = 17.3 - 19.6 mag, which 
includes the brightest red giant of Pal 14 and goes 
down to the limit of faint stars observable with 
UVES. Figure [1] shows Pal 14's color-magnitude 
diagram. 15 of our targeted stars are probable red 
giants and one of the targets may be an AGB or 
evolved horizontal branch star. The significance of 
this different evolutionary state will be discussed 
in Section [3T1 

The spatial distribution of our spectroscopic 
targets is shown in Figure [H The targeted stars 
lie mainly within two core radii with two stars in 
the cluster's outer region. 




Fig. 1.— Col or-magnitude diagram of Pal 14 from 
lHilkerl(|2006l) . The observed targets for the radial 
velocity measurements with UVES are drawn as 
filled and open (red) circles. The open circles de- 
note stars that were subsequently found to be the 
non-members (according to their radial velocity). 
The (blue) filled and open triangles are the stars 
observed with HIRES, the open triangles are the 
non-members. The (green) diamonds are the stars 
observed with UVES and HIRES. Dark grey dots 
are stars within 1 half light radius of Pal 14. 



We used the Besangon Galaxy model (jRobin et al 
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20031 ) to estimate the number of foreground stars 
in our sample. We extracted stars towards Pal 14 
in an area on the sky covering 20r/i, where 
rh = 1.28' is the half-light radius of Pal 14 
(H06). The area covered with our sample stars 
is ^ 2r/i. We selected only those stars located 
in the gray curves shown in Figure [T] and hav- 
ing apparent magnitudes V < 20 mag, and colors 
(B-V) > 0.65 mag. The resulting number of fore- 
ground contaminants in the actual area covered 
predicted by the model is 1. 

UVES was used in its RED 580nm setting cov- 
ering the wavelength ranges 476 - 577 nm (in 
the lower chip) and 584 - 648 nm (in the upper 
chip). We divided the 16 target stars into three 
setups according to their brightness: the bright 
setup, containing the five brightest stars in the 
magnitude range my — 17.37 - 18.52, was ob- 
served for 4 X 60 min in total. The medium setup 
with the four next fainter stars (my = 18.56 - 
19.05) was observed for 6 x 60 min in total The 
faint setup , which included the seven least lumi- 
nous stars (my = 19.19 - 19.76), was observed for 
11 X 60 min in total. The observations were carried 
out in service mode within two observation peri- 
ods, between May 30, 2006, and March 27, 2007. 
The pipeline reduced spectra (R = 60 000) were 
used for the subsequent analysis. 

The zero points in the reduced spectra were not 
identical. The sky emission lines in the single Ih- 
exposures were shifted with respect to each other. 
To correct for this we shifted the spectra to a com- 
mon position of the sky emission lines. As a sky 
zero point location we used the sky lines in one of 
our own observed sky spectra, which we defined as 
reference spectrum. The resulting, shifted science 
exposures were further corrected for the heliocen- 
tric velocity shift. Finally all the shifted single Ih- 
exposure spectra were co-added for each star. In 
this way we get for the brightest star a S/N ~ 16 
and for the faintest star S/N = 4. 

2.2. Spectroscopy v^fith HIRES 

Within a program to study the internal kine- 
matics o f outer halo GC s (for details of the pro- 



gram see Cote et al. I2OO2I) spectra for 11 candidate 



red giants in the direction of Pal 14 were obtained 
using the High Resolution Echelle Spectrograph 
(HIRES; Vogt et al. 1994) mounted on the Keck 
I telescope. The spectra, which were collected 
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Fig. 2. — Spatial distribution of the spectroscopi- 
cally observed Pal 14 stars. The observed targets 
for the radial velocity measurements with UVES 
are drawn as (red) filled and open circles. The two 
open squares mark the non-members (according to 
their radial velocity). The (blue) open and filled 
triangles are the stars observed with HIRES. The 
open triangles are the non-members. The (green) 
diamonds are the stars observed with UVES and 
HIRES. Dark grey dots are probable member stars 
according to their position in the CMD (see Fig- 
ure [T]). The grey concentric, solid circles are from 
inside out the core radius, the half-light radius and 
the tidal radius (H06); the dashed circle is the tidal 
radius calculated by BGK05. The (magenta) area 
is the HST/WFPC2 coverage. 
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during a single night in May 1998, have a res- 
olution of R = 45 000 (for the 0.866" entrance 
slit) and cover the wavelength range between 
506 - 530 nm. The pro gram stars were selected 
from C MDs published by Harris & van den Bergh 
(|1984) and [Holland fc Halr is (199j). The expo- 
sure times were adjusted on a star-to-star basis, 
varying between 900s and 2400s with a median 
value of 1800s. The spectra were reduced entirely 
within the IRAF0 enviro nment, in a manne r iden- 



tical to that described in Cote et al. I (|2002h 



2.3. Radial velocity 

To measure the radial velocity of our targeted 
stars we cross-correlated our final UVES and 
HIRES spectra with two high-resolution spectra 
of the UVES Paranal Obse rvatory Project (UVES 
POP: iBagnulo et ahlliool : HD37811 (a G7 red 
giant) and HD45415 (a G9 red giant). The cross- 
correlation was done with the IRAF task fxcor. 
The heliocentric radial velocities of our two stan- 
dard stars are: wj/dstsu = (-4.68±0.11) km s~^, 
VHDi5ii5 = (52.70±0.04) km s^^ (Melo 2007, pri- 
vate communication). 

We determined the velocity shift of our sample 
stars relative to each of the two UVES POP stars. 
The UVES camera consists of two CCDs. For each 
pair of a UVES science target star and of a UVES 
POP star, we determined two radial velocities, one 
for the upper UVES chip and one for the lower 
UVES chip. These two velocities are averaged to 
a final velocity relative to the UVES POP star. 
Comparing the relative velocities measured for the 
two UVES POP stars, we find a mean difference 
of 0.1 km s^"'^. Within the errors the two veloc- 
ities are equal. The UVES science stars' radial 
velocity is the mean of the t wo velocities weighted 
by the Tonry-Davis R value (|Tonrv fc Davislll979f ) 
determined by fxcor. 

For the HIRES sample, we determined the ve- 
locity of each science star relative to both UVES 
POP stars. Comparing the two relative velocities 
we find a mean difference of 0.07 km s~^. The 
HIRES stars' radial velocity is the mean of these 
two measured velocities weighted by the Tonry- 



Davis R value. 

For 6 stars we have both UVES and HIRES 
spectra. To determine a common zero point of the 
two different samples we compared the measured 
velocities for these 6 stars. A mean velocity shift of 
Av = 0.64 km s^^ was found. The shift is proba- 
bly due to a different instrumental zero point. The 
final HIRES velocities are corrected for this shift. 
The shift can also be due to binarity or stellar 
variability. For two of the six stars we also have 
UVES measurements at two epochs, within the 
errors the velocities agree very well. Short-period 
binarity and variability can be excluded for these 
two stars. The error of the HIRES measurements 
for the five fainter stars is comparable to the mean 
shift. Five stars have a positive velocity shift and 
only one a negative. If all stars were binaries we 
would not expect a clear spread around a positive 
shift. 

The final radial velocity for the 6 stars, with 
UVES and HIRES spectra, is the weighted mean 
of the measured velocities. For the remaining 15 
stars we only have measurements of one instru- 
ment, therefore this velocity is taken as the final 
radial velocity of the star. 

2.4. Photometry 

We used imaging data obtained with the Hub- 
ble Space Telescope/ Wide Field Planetary Cam- 
era 2 (HST/WFPC2) from the HST archive to 
obtain a deep CMD of Pal 14. The data were 
obtained as part of the proposal GO-6512 (PI: 
Hesser). The same data were used bv lDotter et al. 
(|2008 t). The WFPC2 images cover the entire 



^IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 



area within the cluster's core radius (H06), about 
67% of the area within the nominal half-light ra- 
dius (H06), and only 7% of the area within the 
tidal radius (H06) (see Figure [2]). The pipeline- 
reduced FITS files were run thro ugh multidriz- 
zle/tweakshifts (Koekemoer et all 12002) to refine 
the image registration. All further processing was 
done on the original files together with the re- 
fined shifts, us ing the W FPC2 photometry pack- 
age HSTphot (jDolphin I [2000) and following the 
strategy outlined in the HSTphot User's Guide 
for preprocessing, photometry and artificial star 
tests. As in each subset of well aligned images 
in the same filter, the exposure times differed 
significantly, no co-adding was done. In Fig- 
ure [U we show the CMD of all stars brighter than 
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28*'' mag detected by HSTphot with the foUow- 
ing selections: the HSTphot sharpness parame- 
ter {\sharpness\ < 0.2), HSTphot type parameter 
{type < 3, i.e. the star is either a single star or a 
possible unresolved binary), and magnitude errors 
(T,nag < 0.2 mag. 

To determine the photometric errors we in- 
serted artificial stars with known magnitudes. The 
deviations of the subsequently measured magni- 
tudes to the inserted values let us determine the 
photometric errors shown in Figure [3] 

2.5. Completeness 

For a detailed analysis of the stars in Pal 14, 
we performed radius-dependent artificial star tests 
within HSTphot to determine the completeness of 
the observations. For the artificial star experiment 
we added ~160000 stars onto the image. For 7 
annuli of a width of 0.3 arcmin we counted the 
number of artificial stars retrieved from the image 
with a magnitude not more than 0.2 mag different 
from the input value. In Figure |4] we show the 
seven completeness profiles (gray curves), which 
essentially fall on top of each other. Therefore, 
no radial dependence is observed, which is mainly 
due to the low density of Pal 14. The profile of the 
outermost annulus (solid line with squares) shows 
a decline at slightly brighter magnitudes. This is 
an artificial effect. The number of stars in this an- 
nulus is only 10% of the average number of stars in 
the other annulii. We used an averaged complete- 
ness profile in our analysis, shown as the black line. 
The 50% completeness limit is reached at 771555 — 
27.21 mag. 

3. Spectroscopic Results 

3.1. Individual stellar radial velocities 

In Table [U we list the measured heliocentric 
radial velocities {v(UVES) and v(HIRES)) and 
their errors {crv(uvES) and (Jy{HiRES)) for the 21 
stars in our sample. The listed velocities Vrad are 
the weighted mean of the UVES and the HIRES 
observations. The listed errors are the propagated 
errors, weighted by the Tonry-Davis R value from 
the cross-correlations. Star 4, Star 10, HV051 and 
HV086 all have significantly different velocities 
than the majority of the measured stars which 
are centered around ^ 72.2 km s^^: vstar4 = 
(-32.14 ± 0.16) km s-\ vstarw = (50.44 ± 




Fig. 3. — Color-magnitude diagram of Pal 14. We 
show the remaining stars after applying selections 
in the HSTphot parameters: sharpness, magni- 
tude errors, and type. The CMD contains 2752 
stars. 
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Fig. 4. — Completeness for the two filters F555W 
(lower plot) and F814W (upper plot). For each 
filter the completeness profile for seven annuli of 
width 0.3 arcmin are plotted as gray lines with dif- 
ferent symbols, the black line is the overall com- 
pleteness used in the data analysis. 
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Table 1 

Heliocentric radial velocities of our sample stars 





Star*" 


q(2000) 


(5(2000) 


my 


B - V 


v(UVES) 




V (HIRES) 


f{H 1 RES) 


•"Tad 


rad 


m? 






hh:mm:ss.ss 


O.J.I? 


mag 


mag 


km s^^ 


km s^^ 


km s^^ 


km s^^ 


km s^^ 


km s 




1 




16:11:05.81 


14:57:45.1 


17.37 


1.39 


72.53 


0.07 






72.53 


0.07 


Y 


2 


HV025 


16:10:58.73 


14:56:48.7 


17.77 


1.29 


72.76 


0.09 


71.49 


0.30 


72.47 


0.14 


Y 


3 




16:10:54.90 


14:58:36.7 


18.23 


1.07 


71.75 


0.14 






71.75 


0.14 


Y 




HV051 


16:10:59.98 


14:58:30.1 


18.23 


1.48 






-73.77 


1.17 






N 


4 




16:11:04.98 


14:53:32.3 


18.48 


1.10 


-32.14 


0.16 










N 


5 


HV007 


16:10:59.24 


14:57:22.5 


18.52 


0.99 


71.68 


0.18 


73.23 


0.53 


72.21 


0.30 


Y 


6 


HH244 


16:10:53.36 


14:56:45.4 


18.56 


1.04 


72.58 


0.18 


72.79 


0.46 


72.65 


0.27 


Y 


7 


HH201 


16:10:54.04 


14:57:05.6 


18.70 


1.03 


72.62 


0.18 


72.68 


0.49 


72.64 


0.27 


Y 


8 


HV043 


16:10:56.90 


14:57:56.5 


18.84 


0.99 


71.56 


0.21 


70.97 


0.47 


71.38 


0.31 


Y 




HV086 


16:11:02.66 


14:56:41.1 


18.84 


1.21 






-155.31 


0.85 






N 




HV055 


16:10:58.31 


14:58:26.2 


19.02 


0.86 






73.62 


0.89 


73.62 


0.89 


Y 


9 


HV104 


16:11:01.40 


14:57:60.0 


19.05 


0.99 


73.49 


0.21 


73.53 


0.91 


73.50 


0.43 


Y 




HH042 


16:10:53.20 


14:58:12.0 


19.09 


0.98 






71.94 


0.35 


71.94 


0.35 


Y 




HV004 


16:10:58.03 


14:57:25.1 


19.16 


0.80 






73.23 


0.56 


73.23 


0.56 


Y 


10 




16:11:05.89 


14:58:43.2 


19.19 


0.93 


50.44 


0.19 










N 


12 


HV074 


16:10:56.21 


14:56:32.7 


19.41 


0.90 


71.83 


0.23 






71.83 


0.23 


Y 


13 


HV075 


16:10:56.98 


14:56:25.8 


19.44 


0.92 


72.33 


0.41 






72.33 


0.41 


Y 


14 


HV006 


16:10:59.24 


14:57:19.7 


19.50 


0.90 


71.80 


0.27 






71.80 


0.27 


Y 


15 


HV042 


16:10:55.84 


14:57:43.4 


19.60 


0.69 


69.99 


0.38 






69.99 


0.38 


Y? 


16 




16:10:59.62 


15:01:32.9 


19.68 


0.81 


72.14 


0.43 






72.14 


0.43 


Y 


17 


HV021 


16:11:00.58 


14:56:59.1 


19.76 


0.86 


72.39 


0.32 






72.39 


0.32 


Y 



■^Identification from IHillierl I I2006I) 

'='Idontification from [Harris fc van den Berghl II1984I') and IHolland fc HarrTi l ll992f) 



0.18) km s-\ VHVQ51 = (-74.41 ± 1.17) km 
and VHV0&6 = (-155.95 ± 0.85) km s'^. These 
four stars are categorized as non-members (open 
circles and triangles in Figures [1] & [2]) . The re- 
maining 17 stars are considered to be members of 
Pal 14 (see last column in Table [1]). The measured 
velocity of Star 15 is more than Scr away from the 
mean of the other member stars. Therefore, we 
present all our results including and excluding 
Star 15. 

3.2. The systemic velocity and the velocity 
dispersion 

First, we determined the mean velocity and 
the global velocity dispersion for the two different 
measurement sets, respectively. We used the max- 
i mizat ion method described in IPrvor fc Mevlan 
(|l993l) . The mean velocity for the HIRES mea- 
surements is (72.46 ± 0.29) km s~^ and the ve- 
locity dispersion (0.66 ± 0.26) km s^^. For the 
UVES measurements we find a mean velocity of 
(72.28 ± 0.15) km s~^ and a velocity dispersion 
of (0.50 ± 0.11) km s^i if we exclude Star 15. 
Including Star 15, we find (72.12 ± 0.20) km s^^ 
and (0.70 ± 0.15) km s~^. The measurements of 
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Fig. 5. — Radial distribution of stars with veloc- 
ity measurements in Table [TJ Black crosses indi- 
cate the 17 member stars of our sample of Pal 14. 
Star 15, for which the measured velocity is suspi- 
cious, is labeled. The horizontal solid line marks 
Pal 14's global radial velocity without Star 15, 
and the dotted line the radial velocity including 
Star 15. Further the core and half-light radii are 
indicated by dashed grey lines for an easier com- 
parison with Figure 1111 
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Table 2 

Radial velocity and velocity dispersion of Pal 14 



Instrument 


velocity 


error^ 


dispersion 


errord 




km s^^ 


km 


km s~^ 


km s~^ 


UVES" 


72.28 (72.12) 


0.15 (0.20) 


0.50 (0.70) 


0.11 (0.15) 


HIRES" 


72.46 


0.29 


0.66 


0.26 


combined 


72.28 (72.19) 


0.12 (0.18) 


0.38 (0.64) 


0.12 (0.15) 



'^The first value is without the measurement of Star 15. The value in 
parentheses ineludes the measurement of Star 15. 



the two samples agree very well. 

Second, to determine the overall mean veloc- 
ity and the global dispersion for all stars we also 
used t he maximization method of lPrvor fc Mevlan 
(|l993l ). Including Star 15, we measured a 
mean heliocentric radial velocity for Pal 14 of 
(72.19 ± 0.18) km s'^, excluding Star 15, 
(72.28 ± 0.12) km s"^. Within the error bars 
the two values agree. Our results confirm the ear- 
lier measurements by Armandroff et al. ( 19921 ). 

Figure [5] shows the radial profile of our mea- 
sured velocities (Star 15 is labeled). The cluster's 
mean velocity (for both cases) is marked by the 
solid (without Star 15) and dotted (with Star 15) 
horizontal line. In Table [2] we summarize the ra- 
dial velocity and velocity dispersion measurements 
for the two instruments and for the combined stel- 
lar sample. 

The global linc-of-sight velocity dispersion for 
Pal 14 with Star 15 included is (0.64 ± 0.15) km s'^ 
with 99% confidence limits of 0.41 km s~^ and 
1.10 km s^^. Without Star 15, the hnc-of-sight 
velocity dispersion is (0.38 ± 0.12) km s^^ 
with 99% confidence limits of 0.26 km s^^ and 
0.67 km s-^ Within the errors the two val- 
ues would agree. The theoretical prediction for 
the velocity dispersion of BGK05, for which a 
M/L = 2 was assumed, is ctmond — 

1.27 km s-i 

and CTNewton = 0.52 km s~^. For both cases, when 
Star 15 is included or excluded, our results are 
more consistent with the classical Newtonian pre- 
diction, while the MONDian prediction is outside 
the 99% confidence limits. 

As described above the measured velocity of 
Star 15 seems to be deviant. There are several 
possible explanations for this discrepant velocity 
of Star 15: i) Star 15 is a normal member of Pal 14. 
We performed a Monte Carlo simulation in order 



to evaluate how likely the measured radial veloc- 
ity profile is. In the Monte Carlo simulation we 
randomly drew velocities from a Gaussian distri- 
bution, which was newly initialized for each draw 
by calculating the mean velocity and standard de- 
viation of our measured radial velocities randomly 
convolved with their errors. The radial distribu- 
tions of all draws were added. We performed a 
KS-test of the simulated velocity distribution with 
the distribution of the actually measured veloc- 
ity. The KS-test revealed a < 1% probability that 
the distribution that includes Star 15 comes from 
a Gaussian distribution, whereas the probability 
was ^ 50% that the distribution without Star 15 
is Gaussian. This argues against Star 15's mem- 
bership in Pal 14. ii) Star 15 is not a red giant, 
but more likely an evolved horizontal branch (HB) 
star or an AGB star judging from its position in 
the CMD (see Figure [ij . The used templates of 
a G7 (HD37811) and a G9 (HD45415) red giant 
may not be appropriate for Star 15. in) Star 15 
could be a binary. For our faint UVES sample 
(Stars 10-17) we have observations at two epochs: 
June, 2006 and March, 2007. Theoretically this 
allows us to measure a possible change in veloc- 
ity due to binarity. The faintness and the there- 
fore low S/N of Star 15's spectra does not allow 
us to accurately measure the individual radial ve- 
locity for both epochs. The two measured ve- 
locities are W2006 = (70.64 ± 0.63) km s~^ and 
W2007 = (69.13 ± 0.75) km a~'^. Within the er- 
rors the two velocities are the same. Nonethe- 
less, this does not allow us to exclude long-period 
binarity. iv) A further cause for the large off- 
set of Star 15's velocity could be strong atmo- 
spheric variability, which can occur among AGB 
stars. However, from its position in the CMD, 
Star 15 would be an early-AGB star. In this 
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early pha se, AGB stars are not yet pulsating very 
strongly (jHabing fc OlofssonlEioi ). With essen- 
tially only one observing epoch it is impossible to 
know about the star's variability, v) Another op- 
tion might be that Star 15 is not a member of 
Pal 14. We computed a model velocity distribu- 
tion of stars which are located within the light 
gray curves sho wn in Figure [D usi ng the Besangon 
Galaxy model (jRobin et al ][2qq3) as described in 
Section 12.11 The expected velocity distribution, 
for stars with radial velocities > —160 km s~^, is 
shown in Figure [H The number of stars in each 
bin is scaled to an area of ~ 2r/i, in order to repro- 
duce the actually observed area. We expect about 
9 stars to fulfill the photometric constraints. 0.5 
stars have a radial velocity between 50 km s^^ and 
75 km s~^. Therefore, Star 15 could be a fore- 
ground contaminant. 

4. Photometric results 

In order to make predictions for the velocity 
dispersion in Newtonian and MONDian dynam- 
ics, we first have to determine the mass of Pal 14. 
The measured low velocity dispersion is in excel- 
lent agreement with the theoretical prediction of 
classical Newtonian dynamics, and a very strong 
indicator against MOND. The theoretical calcula- 
tions by HBK09 show the dependence of velocity 
dispersion and mass (see Figure 8 in HBK09) for 
classical dynamics and MOND. For a given veloc- 
ity dispersion the necessary mass is always smaller 
in MOND than in classical dynamics. Our derived 
low velocity dispersion is explainable in MOND if 
we find a low total mass for Pal 14. To constrain 
the mass in Pal 14, we analyzed Pal 14's CMD and 
main-sequence mass function. 

4.1. Color-magnitude diagram 

Figure El shows the HST CMD of Pal 14 with 
the remaining stars after the HSTphot parame- 
ter cuts (see Section for details). The CMD 
reaches ^4 mag below the main-sequence turnoff, 
"co-MSTO = 23.63 ± 0.01 mag, which allows us to 
theoretically determine the cluster's mass function 
down to ~0.49 Mq (see Section SH). The CMD 
shows a well-populated main-sequence (MS), sub- 
giant branch, red giant branch, red horizontal 
branch (HB), and some probable blue straggler 
stars. The presence of a red HB and its impli- 
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Fig. 6. — Expected velocity distribution based on 
the Besangon Galaxy model. Only stars that lie 
within the grey area in Figure[T]and that have a ra- 
dial velocity larger than —160 km s~^ are counted 
(see text for more details). The dash-dotted, ver- 
tical line marks the systemic velocity of Pal 14, 
~ 72.2 km s-^ 
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Fig. 7. — Color-magnitude diagram of our fi- 
nal sample of Pal 14 stars from WFPC2. The 
light gray line shows the derived cluster ridgeline, 
the dark gray line is the best fitted a-enhanced, 
[a/Fe]=-|-0.2, Dartmouth isochrone with an age 
of 11.5 Gyr, E(m555-mgi4) = 0.063, and (m — 
M)555 = 19.45 mag. 
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cations were discu ssed in ISaraiedinil (|l997l ) and 



Potter et al 



( 2008h . As expected there is only lit- 
tle field star contamination of Milky Way stars 
due to the moderately high Galactic latitude of 
Pal 14 and due to the small field of view of 
WFPC 2. Jud ginR from the TRILEGAL Galaxy 
Model ( Girardi et al.ll2005h . the number of con- 
taminating foreground stars on our WFPC2 im- 
age in the CMD-area covered by Pal 14 is ~ 2. 
The width of the main-sequence which we observe 
is due to a combination of the photometric errors 
and binary stars. 

For our further analysis of Pal 14, we applied a 
stricter selection of our stellar sample. We deter- 
mined the cluster's fiducial ridgeline (see Figure[71 
the light gray line). The ridgeline reproduces the 
mean location of the stellar distribution in the 
CMD. To derive the clus ter's ridgelin e we ad opted 
the method described in iGlatt et all (|2008a[ l. We 
selected all stars within 2cr of the ridgeline and 
added the blue stragglers and the HB stars for our 
final sample. The 2 500 stars in our final sample 
are plotted in Figure [71 

4.2. Age & Distance 

Pal 14 is known to be younger than typical halo 
GCs (ISaraiedinil I1997I : iHilker II2OO6I : iDotter et al.l 
20081 ) at its metalhcity. We derived Pal 14's 
age via isochrone fitting. We used the Dart- 



mouth isochrones (jPotter et al.ll2007l ). which have 
been shown to reproduce the location of the 
MS, subgiant branch, an d red giant branch very 



well (|Glatt et al.ll2008d ). We adopted the pub- 



lished spectrosc opically dete rmined metalhcity of 
[Fe/H] = -1.50 (|Harrislll996h . Distance and red- 
dening were treated as free parameters. A large 
number of isochrones was fitted using different 
combinations of age, distance, and reddening. We 
selected by trial-and-error the isochrone that best 
matched the above derived ridgeline. 

With an a-enhanced isochrone, [a/Fe] = -1-0.2, 
our best fit yields an age of (11.5±0.5) Gyr, 
a reddening of E(m555 — 171^14) ~ 0.063 (corre- 
sponding to EfB-V) = EfmRsn - TO814)/1-2 = 0.05 
([Holtzm an et al.l Il995l )l. and an extinction cor- 
rectcd dist a nce ni odulus of (m — M)^^^^ = 19.25 mag 
Saraiedinil (|l997l ) stated the age of Pal 14 is 3- 
4 Gyr younger than the age of simil ar halo GCs, 
H06 d erived an age of 10 Gyr and iDotter et al. 
POOSI ) determined an age of 10.5 Gyr via a- 



enhanced isochrone fitting. An a-enhancement 
is found for many of the Milky Way GCs (see. 



e.g. ICarney 1 119961 ). Our new age determination 



reduces the offset to other halo GCs slightly. 

From our CMD and the isochrone fit, we find 
a dereddened distance to Pal 14 of (71 ± 1.3) kpc, 
which places Pal 14 a bit closer to the Sun than 
previously thought. In comparison , H06 de r ived a 
distance to Pal 14 of 74.7 kpc. Dotter et all (|2008f ) 
derived an even larger distance of 79 kpc. 

4.3. Luminosity function 

The cluster's MS luminosity function was de- 
rived by counting the number of stars, fainter than 
the MS turnoff at mMSTO,o = 23.44±0.01 mag, in 
0.5 mag wide bins separated by 0.1 mag along the 
dereddened m^55 axis. 

Furthermore, the WFPC2 images do not cover 
the entire projected spatial extension of the cluster 
on the sky. Our data cover the entire area within 
the cluster's core radius (H06), about 67% of the 
area within the nominal half-light radius (H06), 
and only 7% of the area within the tidal radius 
(H06). The correction for the missing coverage 
within the half-light radius was done as follows. 
We derived the luminosity function for the stars 
within the annulus between the half-light and the 
core radius (nannuius)- Wc then corrected each 
magnitude bin of the entire distribution propor- 
tionally to the distribution of stars within the cov- 
ered annulus: 



N, 



obs 



''annulus \ 



,A 
Ac 



annulus 



1),(3) 



where Acovered is the area of the annulus covered 
by the WFPC2 image, and Aannuius the area of 
the annulus itself. The final correction was done 
for the photometric incompleteness (see Figure |4] 
and Section 12. 5p . We did not correct for poten- 
tial foreground contaminants. The TRILEGAL 
Galaxy model only predicts a very small number 
of stars on our main-sequence. In Figure |S] we 
show the resulting luminosity function of Pal 14. 
The solid line shows the final number of stars per 
0.5 mag bin. The errors are given as N^^^. 



Dotter et al.l (j2008r ) report an unusual flat lumi- 
nosity function for Pal 14 between V =23 mag and 
28 mag. Their data was not corrected for incom- 
pleteness. Our MS luminosity function shows the 
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same flat behavior, correcting for incompleteness 
does not change the slope dramatically. 

4.4. Mass function 

The function dN/dra oc m~" describes the 
number of stars in the mass interval [m, m + dm]. 
We obtained such a mass function for Pal 14's MS. 
The upper boundary of the MS is at the turnofF, 
niMSTO.o — 23.44±0.01 mag. Usin g the mass es 



function slope range of a = —0.9 ... 1.3 (see Fig- 
ure 1 in iDe Marchi et all (|2007t )). The derived 



given by the 11.5 Gyr isochrone bv iDotter et al 



(120071 ). we have stellar masses on the MS covered 
by our photometry between 0.49 Mq and 0.79 Mq. 
We binned the masses linearly into 10 bins of equal 
width of 0.03 Mq. In Table[21 we list the center of 
the mass bins in the first column, and the number 
of observed stars (Nobs) for each bin in the second 
column. 

We corrected the number of stars per mass bin 
for the same effects as in the case for the lumi- 
nosity function. First, the observed number of 
stars per mass bin was corrected for the miss- 
ing area coverage in the same way as described 
above (Table [3l column 3). Second, the mean 
of the stars' incompleteness was used as a cor- 
rection factor. The corrected number of stars 
per mass bin is listed in Table H column 4. To 
fit a slope to our data we only considered data 
with a completeness factor > 0.50 (see last col- 
umn in Table [3]). This restriction leads to a MS 
mass function covering the range from 0.525 Mq 
to 0.79 Mq. We fitted a slope to our data points 
in log{number) vs log{mass) space. In Figure [9] 
we plot the resulting mass function and the fit- 
ted slope of a = 1.27 ± 0.44 as the gray line 



Potter et al.l (j2008[ ) find a similar mass function 
slope of a « 1.2. The canonical Kroupa IMF 



(|Kroupall200lfl in this notation is 2.35 for the given 
mass range. In Figure [H the observed mass func- 
tion is shown as the dash-dotted line. The dotted 
line is the mass function after the area corrections. 
The solid line denotes the completeness corrected 
number of stars. 

Compared to the canonical slope of 2.35 
Pal 14 's mass function i s flatte r in the given mass 
range. iDe Marchi et al.l (|2007l ) compiled the mass 
function slope in the stellar mass range 0.3 to 
0.8 Mq for 20 Galactic GCs of different sizes, con- 
centrations, positions in the Galaxy, etc. Pal 14 
has a (measured) concentration of c = 0.85 (HOG). 
Clusters with a similar concentration span a mass 



slope is comparable with the slope of similar clus- 
ters. E.g., NGC6809 has a concentration of 0.76 
and a mass function slope of 1.3. This slope 
was derived around the cluster's half-light radius, 
where the impact of mass se gregation is negligible 
(|Paresce fc De MarchilboOOl ). For Pal 14 we see 
an increasing number of stars per unit mass down 
to 0.525 Mq. a sudden decrease below this low- 
mass limit would be a unique case as no Galactic 
GC is known to show an initial rise followed by a 
decrease. 

In principle there are two reasons for such a de- 
pleted mass function: Either Pal 14 did form with 
only few low-mass stars, or the cluster is mass seg- 
regated and lost most of its low-mass stars through 
interaction with the Galactic tidal field. The small 
area covered by the WFPC2 image does not allow 
us to estimate the amount of mass segregation. 
In an upcoming paper wc will discuss the issue of 
mass segregation in Pal 14 based on imaging data 
we obtained at the VLT. 

4.5. Total mass & mass-to-light ratio 

To estimate the mass of Pal 14 we corrected 
for the missing area within the half-light radius. 
We measured an observed mass for Pal 14's main- 
sequence Mms,obs = (1 340 ± 50) Mq (above the 
50% completeness limit). The errors are prop- 
agated from the measured photometry. Tak- 
ing into account the stars brighter than the MS 
turn off, correcting for the missing area within 
the half-light radius and the completeness we get 
Mcor- = (2 200 ± 90) Mq within the mass range 
0.525 Mq to 0.83 Mq. If we extrapolate by as- 
suming that the measured slope of a = 1.27 holds 
down to 0.5 Mq and assume a Kroupa-like mass 
function, a = 1.3 for masses between 0.1 Mq and 
0.5 Mq, we have a total mass within the half-light 
radius for Pal 14 of Mtotjd = (6 020 ± 500) Mq. 

The slope of the mass function for st ars with 
mass e s < 0.5 Mp, is s till under debate ( Kroupal 
2OO2I : lElmegreenl l2008l ). Pal 14 is very far from 
the Milky Way. It may have an eccentric or- 
bit that would bring it much closer to the Milky 
Way at perigalacticon possibly leading to strong 
tidal interaction a nd to an en hanced loss of very 
low-mass stars. Ri cher et al.l (2004, 2008) stud- 
ied the main-sequence mass function of the GCs 
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Table 3 
Mass function of Pal 14 



Bin center 
M@ 




^corrected 


N/ 




Completeness 


0.51 


114 


152 


706 


59 


0.21 


0.54 


249 


330 


542 


30 


0.61 


0.57 


196 


255 


316 


20 


0.81 


0.60 


258 


328 


372 


21 


0.91 


0.63 


201 


274 


282 


17 


0.97 


0.66 


226 


306 


311 


18 


0.98 


0.69 


219 


282 


286 


17 


0.99 


0.72 


212 


271 


273 


17 


0.99 


0.75 


213 


268 


270 


17 


0.99 


0.78 


225 


301 


302 


18 


1.00 


0.81 


289 


366 


367 


20 


1.00 



Note. — Column 1 lists the center of our mass bins, column 2 the num- 
ber of observed stars per bin, column 3 the number of stars per bin after 
correcting for the missing area coverage, column 4 contains the final 
number of stars per bin after correcting for completeness, column 6 lists 
the propagated error on the final number of stars per bin, and column 7 
lists the average completeness value for the mass bin. (The numbers in 
columns 3,4, and 5 are rounded to the nearest integer.) 




Fig. 8. — Luminosity function of Pal 14's MS. 
The dashed line is the number of observed stars 
corrected for the missing area coverage The solid 
line is the number of stars after the correction for 
photometric incompleteness (the grey dots mark 
points with a completeness <50%). The horizon- 
tal bars are the N^^'^ errors. 



NGC 6397 and M 4 down to the hydrogen-burning 
limit. In the cluster cores they found mass func- 
tion slopes of a = —0.7, these cluster centers lack 
low-mass stars. Therefore, we also calculated the 
mass in Pal 14 for a mass function with a linearly 
declining slope for masses < 0.5 Mq towards less 
massive stars {a = —1.0). In that case the lower 
limit for the total mass within the half-light radius 
of Pal 14 is Mtot,hi = (2 930 ± 130) Mq. 

If we assume that light traces mass, then the 
half-light radius will also be the half-mass radius. 
Therefore we double the above numbers to esti- 
mate the total mass of Pal 14. The extrapolation 
with a Kroupa-likc IMF for stellar masses between 
0.1 Mq and 0.5 Mq yields a total mass of Pal 14 
of Mtot ~ 12040 Mq. With the declining mass 
function for masses < 0.5 Mq, wc get a total mass 
of Mtot ~ 5 860 Mq. Considering stellar remnants 
will increase the mass further. 

Using the total mass of Pal 14, we derive the 
mass-to-light ratio. The extrapolation with the 
Kroupa-like mass function yields M/L = (2.2 ± 0.4) Mq/Lq. 
The extrapolation with the declining mass func- 
tion gives M/L = (1.1 ± 0.1) Mq/Lq. 
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5. Discussion 
5.1. MOND? 

In HBK09, we calculated the global line-of-sight 
velocity dispersion of isolated and non-isolated 
stellar systems in MOND for circular orbits. For 
details on the simulation see HBK09. In FigurefTOl 
we plot the two curves from these calculations 
showing the global line-of-sight velocity dispersion 
as a function of stellar mass for the classical (open 
squares) and the modified Newtonian case (open 
circles). For a given total mass the velocity dis- 
persion in the MONDian case is larger than in the 
classical theory. In our case, we observed a line-of- 
sight velocity dispersion (shown as the horizontal 
lines) and derived the cluster's mass (shown as the 
vertical lines). 

We measured a line-of-sight velocity dispersion 
of (0.38 ± 0.12) km s-\ not including Star 15. 
For such a low dispersion, the theoretically pre- 
dicted mass in MOND is 950j:4oo and in clas- 
sical dynamics 8 2001^^^;^ Mq. We have observed 
a lower hmit of (2 200 ± 90) M© (marked in Fig- 
ure [10] by the vertical line labeled observed) con- 
sidering only the area within the half-light radius 
of Pal 14. Already the lower limit excludes the 
MONDian case, as we have observed more stel- 
lar mass than MOND predicts and the stars out- 
side the cluster's half-light radius are not consid- 
ered yet. The total mass of ~12 500 Mq (in Fig- 
ure [TO] vertical line marked extrapolated (a — 1.3j) 
is several times larger than the MONDian predic- 
tion. A declining extrapolation at the low-mass 
end down to 0.1 Mq gives a total mass of Pal 14 
of 5 860 Mq (in Figure fTOl vertical line marked ex- 
trapolated (a = — Ij), which is also clearly higher 
than the MONDian prediction. The resulting dy- 
namical mass-to-light ratio for the classical New- 
tonian case is M/Ldyn = (1.48tJ;7o) ^q/Lq- 

If we include the measured velocity of Star 15, 
we find a line-of-sight velocity dispersion of 
(0.64 ± 0.15) km s~^. According to the theo- 
retical calculation of HBK09 the cluster mass in 
MOND would be 2 600j:i2oo and in classical 



mass [M 



dynamics 24 000_ 



+ 11000 
10000 



Mq. In this case, the ex- 



trapolated mass is still larger than the predicted 
mass in MOND. Also, for the declining mass func- 
tion for masses < 0.5 A/© the total mass is larger 
than the MONDian prediction. 




-0.25 



-0.2 
log(mass [M ]) 



-0.15 



Fig. 9. — Mass function of Pal 14. The lowest 
(dash-dotted) line is the observed mass function 
for stars with masses between 0.49 M© and 0.80 
Mq. The dotted line shows the mass function af- 
ter correcting the number of stars per bin for the 
missing area coverage of the WFPC2 data. The 
top line also includes the correction for the pho- 
tometric completeness. In gray the fitted slope 
a = 1.27 ± 0.44 is shown. 
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Fig. 10. — Theoretically predicted velocity disper- 
sion as a function of mass. The two black curves 
are the predictions in MONDian dynamics (open 
circles) and in classical Newtonian dynamics (open 
squares). The observed velocity dispersions (and 
the errors) are drawn as the two horizontal lines, 
the light gray without Star 15, dark gray with 
Star 15. The vertical lines mark the observed 
lower mass limit and the two extrapolated lower 
mass limits. 
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Although the measured low velocity dispersion 
is an indication of whether MOND or classical 
Newtonian dynamics is correct, one can think of 
a scenario in which the cluster would be governed 
by MOND but shows at the same time a velocity 
dispersion consistent with the classically derived 
(low) value, n MOND the gravitational force is ef- 
fectively stronger than in classical dynamics. The 
stars in a GC which resides in the MOND regime 
therefore acquire a higher internal velocity, thus 
leading to a shorter dynamical t ime and a faster 
relaxation time for the cluster ( Ciotti fc BinnevI 
l2004l : IZhaoll2005h . Therefore, aheady after only a 
couple of orbits around the Galaxy, Pal 14 would 
have lost a large fraction of its low-mass stars and 
stellar remnants, leaving the cluster enriched in 
stars around the main-sequence turnoff and on the 
red giant branch. In one of Pal 14's perigalactica 
(if it is on an eccentric orbit), the cluster would be- 
come partially unbound and would expand, while 
it still resides in the classical Newtonian environ- 
ment close to the Galaxy. The unbound cluster, 
then, would move further outward on its orbit and 
would eventually drift into the MONDian regime 
in the Galaxy's outskirts. As MOND is 'stickier' 
than classical Newtonian dynamics, the stars are 
bound more strongly again. As a consequence, an 
observer may measure a low velocity dispersion, 
similar to the value derived in classical dynamics. 
At the same time, the cluster's mass is small. For 
such a scenario to be valid, the cluster would have 
to be strongly effected by tidal forces, but should 
not move too far in to be completely destroyed. 
Detailed simulations on the influence of radial or- 
bits on the velocity dispersion in MOND are neces- 
sary. Unfortunately no proper motion is available 
for Pal 14 in order to make any constraints on its 
orbit. 

MOND is not the only modification of clas- 
sical Newtonian dynamics. One other possible 
theory is mo dified gravity (MOG; iMoffatI l2005t 
iMoffat fc Tolh 2008.)- MOG explains/predicts 
galaxy rotation curves, galaxy cluster masses, etc. 
and at the same time produces predictions consis- 
tent with classical dynamics for smaller systems, 
e.g. GCs. MOG predicts little or no observ- 
able deviation from classical Newtonian gravity 
for GCs with masse s of a few times 10^ M© 
(jMoffat fc TothI I2OO8I ). Our result is consistent 
with the classical prediction and can, therefore. 



neither support nor contradict MOG. 

5.2. Velocity dispersion profile and dark 
matter 

It is widely believed th at globular c lusters con- 
tain no dark matter (e.g., MoorelT996h . Their dy- 
namical masses clos ely match the values from pop- 
ulatio n synthesis ( McLaughlin fc van der Marell 

20051 ). The velocity dispersion profile of GCs 

should, therefore, show a Keplerian fall-off. IScarpa et al 
(j2007[ ) studied velocity dispersion profiles of 
six GCs in the Galaxy. For the five high- 
concentration clusters in their study {tuCen, 
NGC 6171, NGC 6341, NGC 7078, NGC 7099) 
they found the predicted fall-off in the inner parts 
of the clusters, but also an unexpected flattening 
in the outer parts. On the other hand, for the 
low-concentration cluster NGC 288 they found a 
more or less flat dispersion profile. For the high- 
concentration clusters, the profiles always flatten 
at a radius where the acceleration is around the 
MONDian limit of oq for a mass-to- light ratio of 1. 
To draw any conclusion about MOND from this 
is rather difficult, as the discussed clusters' total 
accelerations are not below ag and therefore the 
effect of MOND is tiny or even not existent; the 
clusters are all too close to the Galactic center. 

All our stars but one are located within 2.5' of 
the center of Pal 14. We derived the line-of- 
sight velocity dispersion profile with running ra- 
dial bins, each bin containing six stars. Figure flT] 
shows the resulting velocity dispersion profile. Be- 
tween 1' and 1.5' we derived the velocity dispersion 
either including Star 15 or excluding Star 15. The 
lower black squares are the case where Star 15 was 
not included, the upper gray squares the case in- 
cluding Star 15's velocity. For the case excluding 
Star 15, we can see (within the errors) a slightly 
declining velocity dispersion profile. The dashed 
curves in FigurefTTlare the theoretically calculated 
profiles of HBK09. If we compare our disper- 
sion profile to the theoretical predictions we see 
a slow fall-off towards outer radii for both. We 
have observed velocity measurements in the inner 
2.5' ('^SO pc ~ 3.6 core radii). On the other hand, 
Scarpa et all (|2007t ) showed the velocity dispersion 
profile of NGC 288, anot her sparse GC with a con- 
centration of c = 0.96 ( Harris! 19961 ). They de- 
scribe the profile to be fiat out to 4.5 core radii. 
In order to improve the significance of the compar- 
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Fig. 11. — The velocity dispersion profile of Pal 14 
using running bins with six stars in each bin. The 
black squares denote the velocity dispersion with- 
out Star 15. The gray squares denote these bins 
where Star 15 was included. The black and gray 
dashed curves are the theoretical dispersion pro- 
files if Star 15 was included and excluded, respec- 
tively. The vertical, dotted lines are the core and 
half-light radii (H06), respectively. The arrow at 
the top of the plot marks the radial distance of 
Star 15 from the cluster center. 



ison of the theoretical prediction and the observa- 
tional data for Pal 14 as well as of the compari- 
son with similar clusters, spectroscopic data out 
to larger radii are needed for Pal 14. 

We treat this GC the same way as dwar f 
spheroidal (dSph) galaxies in Madau et al. ( 20081 ) 
to calculate the central density, using po = IGGrycr^/r^ MqPC~ 
setting r] = 1, Tc = 0.7' = 14.5 pc, and 
cr = 0.38 ± 0.12 km s'^ We find a central den- 
sity of Po = 0.1 ± 0.07 Mqpc^^. a value which 
is very similar to values found for dSp h galax- 

Madau et all (|200g )). 



les 
On 



[see 
the 



e.g. Table 
other hand, 



1 in 

if we 



derive the density 
within the half-light radius from our mass estimate 
Mtot.M = (6 020 ± 500) Mq and = 1.28' = 
26 pc, we find p = 0.08 ± 0.01 Mqpc-^. Within 
the errors the two values agree. We do not need 
to assume DM for Pal 14. 

6. Summary 

Modified Newtonian dynamics has proven to be 
quite succ essful on galactic and also on intergalac- 
tic scales (|Sanders fc McGaughir2002f) . However, 
not only galaxy size objects must be correctly ex- 
plained by MOND. Objects with similar low accel- 
erations, for which there is no need for additional, 
unseen matter such as GCs must be described cor- 
rectly by this modified theory, as well. Hence, 
we have studied the outer halo GC Pal 14 to test 
whether modified or classical Newtonian dynamics 
applies. Pal 14 has an internal and external accel- 
eration that are both significantly smaller than ag . 
Also, the total acceleration of stars in Pal 14 is still 
significantly smaller than aq and therefore. Pal 14 
is an excellent test object for the two theories. 

We determined the radial velocities of 17 gi- 
ant stars in Pal 14. Using the measurements of 
all 17 giants, we confirmed the cluster's mean 
radial velocity of (72.19 ± 0.18) km s~^ and mea- 
sured a global line-of-sight velocity dispersion of 
(0.64 ± 0.15) km s'^ (see Section E] for details). 
Excluding Star 15, we find a similar systemic ve- 
locity of Pal 14 of (72.28 ± 0.12) km s'^ and a 
lower velocity dispersion of (0.38 ± 0.12) km s^^. 
These velocity dispersions lead to dynamical 
masses of 950^400 Mq in modified dynamics, and 
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4000 



Mq in classical dynamics for the case 



without Star 15. In the case including Star 15 we 



expect total masses of Pal 14 of 2 600 
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The mass function of Pal 14 has a slope of 
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0.78Mq and is thus flatter than the canonical mass 
function. This is consistent with the cluster being 
formed mass segregated with a normal (canoni- 
cal) IMF but suffering major rn ass loss through 



gas expulsion ([Marks et al.ll2008l ). The HST im 
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